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Abstract: We demonstrate novel features in the behavior of the second and third order 
non-linearity parameters of the curvature perturbation, namely, fjyi and gNL, arising from 
non-linear motion of curvaton field. We investigate two classes of potentials for the curvaton 
- the first has tiny oscillations super-imposed upon the quadratic potential. The second 
is characterized by a single 'feature' separating two quadratic regimes with different mass 
scales. The feature may either be a bump or a fiattening of the potential. In the case of 
the oscillatory potential we find that as the width and height of superimposed oscillations 
increase, both /tvl and qnl deviate strongly from their expected values from a quadratic 
potential, /tvl changes sign from positive to negative as the oscillations in the potential 
become more prominent. Hence, this model can be severely constrained by convincing 
evidence from observations that fiy^ is positive, qnl, on the other hand, acquires very 
large negative values. For the the single feature potential, we find that /^vl and gi\iL exhibit 
oscillatory behavior as a function of the parameter that controls the feature. 
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1. Introduction 

The inflationary paradigm has become an important ingredient of modern cosmology. 
Inflation provides a natural explanation for the production of the first density perturbations 
in the early universe which seeded the formation of the large scale structure (LSS) in 
the distribution of galaxies and the temperature anisotropics in the cosmic microwave 
background radiation (CMBR) Q. However, the precise details of the mechanism for 
generating the primordial curvature perturbation is not fully established. The standard 
mechanism is via the quantum fluctuations of the inflaton field. An alternative scenario 
which frees the inflaton from the job of generating perturbations, besides giving rise to 
inflation, is the curvaton scenario |^ |5|. The curvaton is assumed to be a light scalar 
field which begins evolving at the end of inflation. Its energy density is assumed to be 
subdominant during inflation, but it can share a significant part of the total energy in the 
universe before its decay. The entropy perturbations caused by the curvaton field finally 
get converted into adiabatic perturbations. 

A large number of light scalar fields are expected to be present in any fundamental 
theory that goes beyond the standard model of particle physics. During the inflationary 
era these fields would have had the same amplitude of quantum fluctuations. It is plausible 
that at least some of them played important roles in the early universe, for example, as 
the curvaton field. An important distinguishing property of the curvaton scenario as the 
generating mechanism for primordial perturbations, from the standard single slow-rolling 
field picture, is the possibility for the primordial perturbations to have large deviations 
from Gaussian distribution. This property becomes very attractive in the light of the 
recent result from WMAP which suggests that primordial non-Gaussianity may be large 
Q. The non-Gaussianity generated in the curvaton scenario must have a local shape 
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because it is generated on superhorizon scales. Then the curvature perturbation can be 
expanded at the same spatial point to non- linear orders, as, 

C(x) = C,(x) + ^/nl (Cg'Cx) - (Cg'» + ^9NL (C,'(x) - 3(C,')C,) + ... , (1.1) 

where /nl and g^L are the so-called non-Gaussianity parameters. The WMAP 7yr result 
implies a constraint on the size of local form bispectrum as /at^ = 32 it 21 at 1 — a level. 
The limits on gj\[L are —3 x 10^^ < qnl < 8 x 10^^ from LSS and similar limits are 
obtained from CMB data from WMAP Syr data Q as well. A convincing detection of the 
local form non-Gaussianity will rule out all single-field inflation in a model-independent 
way. 

In the simplest case the curvaton potential is assumed to have a quadratic form and 
the typical size of the bispectrum is bounded by the tensor-scalar ratio, as, /nl < lO^r^/^, 
Q. Since the curvaton field evolves linearly in this case, the size of the trispectrum is 
linearly related to that of the bispectrum, as, 

gNL ^ -y/7VL. (1.2) 

There is, however, no reason for the above relation to hold in general from the viewpoint 
of fundamental theory. Apart from the quadratic potential, all the curvaton models that 
have been considered so far in the literature focus on potentials which deviate from the 
quadratic form at large field values but tend to the quadratic form at small field values. The 
predictions of such models, particularly the level of non-Gaussianity, are then compared 
with those from the quadratic potential so as to understand their distinguishing features. 
Clearly, the distinction becomes more prominent as the initial curvaton field value be- 
comes larger and larger. A distinct signature of departure of the curvaton potential from 
quadratic form is a breakdown of the relation (1.2). If the curvaton self-interaction term 
becomes dominant, giving rise to higher order corrections in the curvaton potential, the or- 
der of magnitude of gj^L can be 0{ffj^) |jlO|, 11, 12, |l3[. The predictions of curvaton model 



with nearly quadratic potential are investigated in [0, 15, |T^, 17|, where the non-linear 
evolution of curvaton after infiation but prior to its oscillation is taken into account. An- 
other promising curvaton candidate is the pseudo-Nambu-Goldstone boson - axion, whose 
potential significantly deviates from quadratic form around the top of its potential. A 



numerical analysis of the axion- type curvaton model is discussed in |18, From the 



viewpoint of fundamental theory, one can generically expect multi curvatons models and 
such a model is investigated in ||2^, While the discussion this far has ignored scale 
dependence of the bispectrum and the trispectrum, it is possible that in the future such 
scale dependences may become accessible to experimental observation and hence important 
Other papers of relevance are p^, E3, |2^, E^, |2M] . 



We investigate two new curvaton models different from the ones described above. The 
first is a potential which has tiny oscillations superimposed on the quadratic form. The 
resulting effect on the curvaton evolution is that it experiences the small bumps of the 
oscillations in the potentials during the stage of it evolution when it undergoes oscillations 
about the minimum of the potential. As a consequence, the curvaton evolution during 
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this stage is non-linear (the curvaton equation of motion is not that of a damped simple 
harmonic oscillator), making it significantly different from the curvaton oscillation about 
the minimum of a quadratic potential. Our goal is to calculate the non-linear curvature 
perturbation up to cubic order and obtain the predictions for non-Gaussianity from such 
a model. We find very interesting new implications for the non-linearity parameters f^L 
and giyL arising in this model. First, fjyi is no longer restricted to have positive values. 
Depending on the amplitude and the frequency of the superimposed oscillations on the 
potential, it can take a wide range of both positive and negative, with a switch of sign 
from positive to negative, qnl , on the other hand, remains negative and can take large 
negative values. The sign switch of /nl brings up the possibility that the most important 
contribution to primordial non-Gaussianity could come from the qnl term, with /tvl being 
negligibly small. 

The second model we discuss is a class of potentials characterized by a single feature 
separating two quadratic regimes with different mass scales. The feature depends on a 
single parameter and depending on the sign of the parameter, it can be either a single 
bump or a fiattening of the slope of the quadratic potential at some characteristic scale. 
We find that the effect of the feature on Jnl and g^L is rather dramatic, causing them to 
oscillate with increasing amplitude as the strength of the feature increases. 

This paper is organized as follows: in section 2, we briefly summarize the method 
for computation of the non-linear curvature perturbation using the 5N formalism and the 
curvaton equation of motion. In section 3, we describe the specific forms of the curvaton 
potentials we are considering here and display our results for the non-linear corrections to 
the curvature perturbations. In section 3.1, we discuss the case of the washboard potential, 
while in section 3.2, we discuss the single feature potential. We end with a summary of 
our results and comments in section 4. A brief description of the curvaton with quadratic 
potential is given in the appendix to highlight the differences from our study and novelty 
of our results. 



2. The non-linear curvature perturbation 

On sufficiently large scales, the curvature perturbation on the uniform density slicing can 



be calculated by using the so-called 6N formalism [M, RQ, 31, 32, pS]. Starting from any 



initial fiat slice at time tini, on the uniform density slicing, the curvature perturbation is 

C{t,^) = 6N = N{t,^) - No{t), (2.1) 

where N{t, x) = lna{t, x)/a{tini) describes the local expansion of our universe, and A^o(*) = 
In a (t) / a (tini) is the unperturbed amount of expansion. In curvaton model, the difference 
between local expansion and the unperturbed expansion is caused by the quantum fiuctu- 
ations of curvaton field during inflation. Therefore 

C = N^Ja + ^N.^da^ + ^N^^^^da"^ + ... , (2.2) 
2 D 
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where iV,^ = dN/da, N^„„ = d^N/da'^ and N^„„„ = d^N/da^. Considering 6a = H^/2n, 
the amphtude of the power spectrum generated by curvaton field is 

Pc,=Nl[ll)\ (2.3) 
and the non-Gaussianity parameters are given by 

f^L = (2.4) 

= 54^' ^'-^^ 

where is the Hubble parameter during inflation. On the other hand, the amplitude of 
the tensor perturbation only depends on the inflation scale, namely 

Pt = ^. (2.6) 

Thus the tensor-scalar ratio r is given by 

r^PTlP,, = j^,. (2.7) 

Here wc consider the simplest version of curvaton scenario where the quantum fluctuations 
of convaton field contribute the total curvature perturbation. 
After inflation, the equations of motion are 

= ^iPr + Pa), (2.8) 
Pr + 4:Hpr = 0, (2.9) 

p, = ^a^ + V{a), (2.10) 

a + 3Ha + ^ = 0, (2.11) 
da 

where pr and p^ are the energy densities of radiation and curvaton respectively, and V{a) 
is curvaton potential. In order to numerically solve the above differential equations, we 
define the reduced curvaton field a and reduced curvaton potential V{a) as follows 

a = a/a^, (2.12) 

= (2-13) 

where a* is the vacuum expectation value (VEV) of curvaton field in the inflationary era. 
Now the equations of motion can be simplified to be 



N' 



-4N I ^* f^~'2, \rl~-^ 



1 



(2.14) 



a" + ZN'a' + = 0, (2.15) 

da 
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where N{x) = Ina(t), a = g^^is^a = Hf^-fm'^, and the prime denotes the derivative with 
respect to dimensionless time coordinate x = mt, and the Hubble parameter becomes 

H = mN'. (2.16) 

The solution for the subdominant curvaton with quadratic potential is analytically dis- 
cussed in the appendix. 

The scale factor can be rescaled to satisfy a{tini) = 1, or equivalently N{tini) = 0. 
For numerical calculation, we also need to input the value of a. If the vacuum energy of 
inflaton suddenly decays into radiation, a reasonable choice is a = Hf^j/rm? which is much 
larger than one. However we don't know its value exactly. But as long as a is large enough 
it does not affect our numerical result because the curvaton field almost does not move 
when the Hubble parameter is much larger than its mass. For example, it is reasonable to 
assume that the Hubble parameter at the inflationary era is one order of magnitude larger 
than the curvaton mass at least and then we set a = 10^ in this paper. 



3. The models 



In this section we consider two new curvaton models which have some small features around 
the exactly quadratic form of the curvaton potential. We can expect that these features 
will introduce non-linear effects to the oscillating curvaton field and consequently affect 
the non-Gaussianity parameters. Our aim is to calculate the precise effects. Note that 
these effects are different from what was considered in 16, where the non-linear 



evolution of curvaton after inflation but prior to its oscillation was considered. Since the 
non-linear nature of the curvaton motion makes analytic solutions extremely difficult to 
obtain, we rely on numerical methods to get our results. We solve the Eqs. ( 2.14| ) and 



(2.15) as a coupled set of differential equations for each potential under consideration. 



3.1 Washboard curvaton model 

Let us consider a curvaton potential which has tiny oscillations superimposed on the exactly 
quadratic form. We call it the washboard model and it takes the following explicit form 

V[a) = imV2 + Vo{l- cos(^)) , (3.1) 

where Vq <^V* = ^m^ci^. The reduced potential of a is 

F(a) = V + e(l-cos(<7/5)), (3.2) 

where ^ 

e = ^-3, 5 = —. (3.3) 

Here e measures the size of the correction and 5 characterizes the period of oscillation of 
the correction term in the washboard potential. The reduced potential is shown in the left 
panel of Fig. |l] for easy visualization. When a ^> the potential is almost quadratic. If 
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Figure 1: The washboard curvaton potential given by Eq. (3.2) is shown on the left panel for 
visual comparision with the corresponding quadratic one. The parameter values are e = 5 x 10^'* 
and S = 10~^. We have chosen a large value of e in order to make the oscillations clearly visible. 
The right panel shows the curvaton oscillations about the potential minimum for the quadratic and 
washboard cases, with the same initial field value given by a^/AIp = 0.1. The parameter values for 
this plot are e = 10^^ and S = 10~^ 



<T <C 5, then V{d') ~ ^(1 + e/(5^)a^. Then the curvaton potential is roughly quadratic as 
well, but has a deformed mass. 

The dynamics of curvaton field after inflation is governed by 

a" + Act' + a + I sin {a/5) = 0. (3.4) 
2x 

Even though the correction to the potential is small, the dynamics of curvaton field becomes 
significantly non-linear if the period of the correction term is small enough. Here we 
consider the case in which the dynamics of curvaton is dominated by the mass term in the 
beginning, which implies e/5 < 1. Once the amplitude of the curvaton oscillation drops 
below e/6, the curvaton evolves non-linear ly. On the right panel of Fig. |l| we have plotted 
the oscillation of the curvaton field about the minimum of the potential for the quadratic 
and the washboard potential cases, for the same initial field value given by a-^/Mp = 0.1. 
We can see that the amplitude of oscillation in the washboard case decreases faster than 
the quadratic case. Moreover, the frequency of oscillation for the washboard curvaton is 
time dependent, it oscilates about the constant frequency of the quadratic case. 

We now illustrate how the small features in the curvaton potential play an important 
role for the non-Gaussianity parameters. First we solve for N^^j, N^^^j and N^^^^j and then 
obtain /tvl and gNL from them. In general, we need to scan four independent parameters, 
namely, T/m, a^:/Mp, e and 6 in order to satisfy observational constraints such as amplitude 
of perturbations and the limits on /^vl and gNL- Our strategy here is to fix T/m, a^/Mp 
and 5 and obtain f^L and gNL SiS functions of e. Our results are obtained for two values of 
r/m and 5 each, to understand how these parameters systematically affect Jnl and gNL- 
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Figure 2: N^^y, N^^a, N^rraa, /nl and qnl are shown as functions of e for the washboard potential, 
for fixed values oi 5, T/m and a^jMy. We have shown plots for two different values of F/to in order 
to demonstrate the systematic variation as we change F/m. 



In the case of the quadratic potential T/m is typically required to be of the order of 
10^*^ for the amplitude of perturbations to be COBE normalized. Evolving the equations 
numerically till the energy density of the curvaton decreases to such small value is pro- 
hibitively time consuming. Moreover, for the purpose of capturing the essential features of 
the dependence of /atl and g^q^ on e and (5, it is enough to fix T/m at a relatively larger 
value. To demonstrate this point, in Fig. |2|we have plotted for T/m = 10~^ and 2 x 10~^, 
how A^cr, -/VfTo-, Nfjfjfj, /tvl and qnl vary as functions of e, for fixed values of 5 = 10~^ and 
a^/Mp = 0.1. We can see that T/m systematically changes the amplitudes of Jnl and 
qnLi but does not alter the essential functional shapes. The correctness of the numerical 
calculations are tested by ensuring that in the limit e — 0, Jni and qnl tend to their 
analytically expected values for the quadratic potential, as clearly seen in the figure. As 
e increases, f^L and qnl get strongly affected and deviates from their expectation from 
quadratic potential. f^L crosses over from positive to increasingly negative values as e in- 
creases. On the other hand, qnl remains negative throughout but its magnitude becomes 
very large as e increases. The inset figures in the panels showing A'^o-o-o- and qnl zoom in on 
the e — )• region to show them approaching the negative values expected from quadratic 
potential. 

Next, in Fig. |3| we have plotted N^, N^j^j, Nf^cra, Inl and gNL as functions of e, for 
two different values of 6. We have chosen 6 = 10~^ and 1.8 x 10~^ and fixed T/m = 10~^ 
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and a-t/Mp = 0.1. We see that for very small e, varying 5 has little effect on the behavior 
of fNL and qnl- This can be explained by the fact that e — )• kills off the oscillations 
superimposed on the potential, regardless of the frequency of oscillations which is controlled 
by 5. At relatively larger values of e, the effect of 5 becomes prominent. As the deviation 
of fNL and gNL from the quadratic potential behavior increases as 5 decreases, due to the 
increase in the frequency of the oscillations in the potential. As in Fig. ^ the inset figures in 
the panels showing A'^o-o-o- and qnl zoom in on the e — )• region to show them approaching 
the negative values expected from quadratic potential. 




Figure 3: Same as Fig. 2 but for two different values of (5, with F/to and CT*/Mp kept fixed. 
3.2 Single- feature curvaton model 

In this subsection we consider a curvaton model with the potential given by 

where re > 0, and M is an energy scale which measures the position of the feature. In the 
regime a ^ M or o" ^ M, the curvaton potential has a quadratic form, but around a ~ M 
the potential deviates from quadratic form. As in the previous subsection, we define a 
reduced potential, as follows 
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where 



The reduced potential given by Eq. (|3.6|) is shown in Fig. for n = 2 and d = 0.1. The 



d = ^. (3.7) 




Figure 4: The single feature curvaton potential given by Eq. (3.7) is shown on the left panel for 
comparision with the corresponding quadratic one. d is fixed to be 0.1 and we have plotted for 
c = 2 and —1 to show how the nature of the feature changes with the sign of c. The right panel 
shows the corresponding curvaton oscillations about the potential minimum, with the same initial 
field value given by a^/Mp = 0.1. 

nature of the feature depends on the sign of c. If c is positive, then there is a bump, 
whereas, a negative c changes the slope of the potential to make it flatter around some 
scale set by the parameter d. 

The equation of motion for the reduced curvaton field becomes 

/ / / 1 - (n- l)(^/(f)2"\ ^ ^ 

We restrict our analysis here to n = 2. If d <^ 1 and the initial curvaton field value is large 
enough, then the curvaton evolves linearly prior to its oscillation. We choose d = 0.1. As 
in the washboard curvaton model, we choose Tfj/m to be 10~^ and a^jM^ = 10^^. Then 
we solve for /atl, scanning the parameter c. On the left panel of Fig. ^, we have plotted 
-^,0-) -^,(T(T and /tvl for the single-feature potential. As shown in the figure, f^L oscillates 
about zero with increasing amplitude as |c| increases. We can see Jnl flattening out and 
approaching the expected value from the quadratic potential as |c| — t- 0. On the right 
panel of the same fugure we have plotted N^cr, N^^j^j^j and gNL- We again obtain oscillatory 
behavior of qnl as c varies. Similar to /atl, we can see the curve flattening out near c = 
for gi^L to assumes the value expected from quadratic potential, as shown in the inset 
figure on the right bottom panel. 
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Figure 5: N,a-, N,c,<t, -^.o-o-o-, Inl and g^L are shown as functions of c for the single-feature 
potential, with d, T/m and a^/Mp kept fixed. 



4. Conclusion and discussion 

We have studied two new curvaton models in this paper. The first is the washboard 
model where the potential has tiny oscillations superimposed on the quadratic form, and 
the second one has a potential with two quadratic regimes having different mass scales 
separated by either a bump or a flattening of the potential. For the washboard model we 
have investigated in detail how the two parameters that control the oscillations, namely, the 
amplitude and the frequency, affect the non-linear corrections to the curvature perturbation 
via their effect on f^i and qnl- We have shown that the relation g^i oc —/nl, which 
holds for the quadratic potential, is no longer valid in this case. We also found that there 
is a wide range of both positive and negative values for Jnl, while g^L remains negative 
but its magnitude can be very large depending on the model parameters. In comparision, 
the quadratic potential restricts /nl to be positive and qnl to be negative. For the single- 
feature model we have again calculated /nl and qnl, and demonstrated that they strongly 
depend on the strength of the feature and oscillate as the strength increases. 

What is new in the models considered here is that the curvation motion as it oscillates 
about the potential minimum is non-linear, unlike other models that have been considered 
so far in the literature. The results that we have found have interesting implications 
for searches for non-Gaussianity in observational data. The fact that fi\fL can switch 
sign at some parameter values implies that it is possible that the non-linear contributions 
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to the curvature perturbation could be coming from ^^VL alone, with /nl being close 
to zero. Similar result was obtained in [14| in the context of curvaton potential with 
non- linear corrections to the quadratic term. It is also possible that both f^i and qnl 
contribute comparably with the same or opposite signs. The present work thus throws up 
the need to understand different sources of primordial non-Gaussianity and how they can 
be distinguished in the observational data. It is important to devise observables which can 
distinguish them. Such studies have been initiated in |34, 35 1. We also want to mention that 
f]\[L and gisiL are controlled by two independent geometric quantities which characterize the 
hyper-surface in field space on which multi-field inflation ends. This has the implication 
that all of the possible results for /tvl and g^L in curvaton models can be realized by 
tuning these two independent geometric quantities, as shown in [|^, 37, 

In principle, the three free parameters in the washboard model, T/m,e and 5, can 
be constrained by using the observational constraints on /nl, 9nl and the amplitude 
of perturbations. The parameters c and d in the single-feature model can be similarly 
constrained. However, such a full scan of the parameter space is beyond the scope of 
the present analysis. Our purpose in this paper has been to understand the systematic 
behaviors of /tvl and g^i as functions of the model parameters. We will tackle the problem 
of scanning the parameter space in a future work. 
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A. Curvaton model with quadratic potential 

For the curvaton model with quadratic potential, from p3|, we have 



fNL = 77 - -fn, (A.l) 

4/d 3 6 

25 25 125 „ 25 „2 . 
= "6f^ ^ l08 + + 18^^' ^^-^^ 

where 

After inflation the universe is dominated by radiation and the Hubble parameter is related 
to the cosmic time t by H = ^ . The equation of motion of curvaton field with quadratic 
potential becomes 

a" + ^a' + a = 0, (A.4) 
2x 
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whose solution is 

a = 2i/^r(5/4)x-i/Vi/4(x), (A.5) 

where Ju{x) is the Bessel function of the first kind. Therefore the energy density of curvaton 
is given by 

Pa = {a^ + = ^^mV.^x-V^ (j^^/J,) + ^J,)) . (A.6) 

Adopting the sudden decay approximation, we have 

2 I 

in the hmit ofxD = 5^3>l. In the Hteratures, i^a,D = which is roughly half of 

our exact retult. 
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